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We have studied the interplay between orbital ordering, Jahn- Teller and GdFeOs-type lattice distor- 
tions in perovskite-type transition-metal oxides using model Hartree-Fock calculations. It has been 
found that the covalency between j4-site cations and oxygens causes interaction between the Jahn- 
Teller and GdFeOs-type distortions. The present calculations explain why the d-type Jahn- Teller 
distortion and orbital ordering compatible with it are realized in LaMnOs, YVO3 and YTiOa. 






X3 
O 

o 



>■ 
o- 

o: 

O' 

on; 



O 
o 



X 



I. INTRODUCTION 

Orbital ordering and concomitant Jahn- Teller (JT) 
distortions are observed in some perovskite-type 3d 
transitioitjiujetal compounds such as KCuFa and 
LaMnOa.El Q In the perovskite-type lattice, there are two 
possible JT distortions depending on the stacking of the 
elongated octahedra along the c-axis as shown in Fig. 
n(a).B In the d-type JT distortion, the elongated axes 
of the octahedra are parallel along the c-axis. On the 
other hand, the elongated axes are rotated by 90° along 
the c-axis in the a-type JT distortion. While LaMnOs 
{d'^), YYOf^Lf) and YTiOa (d^) have the d-type JX 
distortion,yili LaVOa (d^) has the a-type JT distortion.Q 
In KCuFa, both the d-type and the a-type JT distortions 
are observed.B Hartree-Fock calculations which consider 
the hybridization between the transition-metal 3d and 
oxygen 2p orbitals predict that the orbital ordered state 
compatible with the a-type JT distortion is lower in en- 
ergy than that with the d-type JT distortion for d^ and 
d^ systems and that the two states are degenerate for 
d* and d^ systems J3 Therefore, one cannot explain why 
the d-type JT distortion is realized in LaMnOa, YVO3 
and YTiOa by considering the energy gain due to orbital 
ordering alone. 

Perovskite-type ABO3 compounds with relatively 
small A-site ions undergo the GdFeOa-type distortion 
which is capsed by tilting of BOe octahedra as shown in 
Fig. 0(b) B While LaMnOs, YVO3, LaVOg and YTiOg 
are accompanied by the GdFeOa-type distortion, KCuFa 
has no GdFeOs-type distortion. In addition, the mag- 
nitude of the GdFeOa-type distortion becomes larger in 
going from LaVOa and LaMnOs to YVO3 and YTiOa. 
Here, one can notice that the compounds with the larger 
GdFeOs-type distortion tend to have the d-type JT dis- 
tortion. It has been pointed out by Goodenough that 
the covalency between the A-site and oxygen ions {A-fl 
covalency) is important in the GdFeOs-type distortion.EHl 
In this paper, we have studied the relationship between 
the GdFeOs-type and JT distortions considering the yl-0 
covalency and explored the reason why orbital ordering 
compatible with the d-type JT distortion are favored in 



LaMnOa, YVO3 and YTiOa in terms of the interaction 
between the two distortions. 



II. METHOD OF CALCULATION 

We have employed lattice models for the perovskite- 
type structure in which the transition-metal 3d, the oxy- 
gen 2p and the A-site cation d orbitals are included. The 
on-site Coulomb interaction between the transition-metal 
3d orbitals, which is essential to make the system insu- 
lating and to cause orbital ordering, are-expressed us- 
ing Kanamori parameters u, u' , j and j'.EJ The charge- 
transfer energy A is defined by f^ — fp + nil, where e^ 
and gp are the energies of the bare transition-metal 3d 
and oxygen 2p orbitals and U {—u-20/9j) is the multiplet 
averaged d-d Coulomb interaction energy. The hybridiza- 
tion between the transition-metal 3d and oxygen 2p or- 
bitals is expressed by Slater-Koster parameters (pdcr) and 
{pdn). The ratio {pda)/{pd-K) is fixed at -2.16.E3 A, U, 
and {pdu) can be deduepd from cluster-model analysis of 
photoemission spectra.El Although the error bars of these 
parameters estimated from photoemission spectra are not 
so small [-- ±1 eV for A and U and ~ ±0.2 eV for {pda)], 
the conclusions obtained in the present calculations are 
not changed if the parameters are varied within the error 
bars. 

In the present model, unoccupied d orbitals of the A- 
site cation such as Y 4d and La 5d are taken into account. 
The hybridization term between the oxygen 2p orbitals 
and A-site cation d orbitals is expressed by {pda)A and 
(pdTr)^. The ratio (pda) a / (pdn) a is also fixed at -2.16. 
The hybridization term between the oxygen 2p orbitals 
is given by (ma) and {ppn) and the ratio (j>pa)/{ppTr) 
is fixed at -4Jl3 It is assumed that the transfer integrals 
ipdcr) and {pdcr)A are proportionsLLto d^'^'^ and (ppcr) is 
to d~^, where d is the bond length.tJ Without the JT and 
GdFeOs-type distortions, the bond length between two 
neighboring oxygens and that between the oxygen and 
the A-site cation are \/2a, where a is the bond length 
between the transition-metal ion and the oxygen, (ppcr) 
and {pda)A for bond length of \/2a are assumed to be 



-0.60 and -1.0 eV, respectively. 

In the GdFeOs-type distortion, the four octahedra in 
the unit cell are rotated by angle of uj around the axes 
in the (0,1,1) plane in terms of the orthorombic unit cell. 
Here, we model the GdFeOs-type distortion by rotating 
the octahedra around the (0,l,0)-axis or the 6-axis [see 
Fig. |l](b)]. The subsequent small rotation around the 
a-axis is required to retain the corner-sharing network 
of the octahedra. The magnitude of the GdFeOs-type 
distortion is expressed by the tilting angle to around the 
&-axis. It is important that the A-site ions are shifted 
approximately along the &-axis to decrease the distance 
from the A-site ion to the three closest oxygen ions and 
increase the distance to the three next closest oxygen ions 
as shown in Fig. |^(b). Here, it is assumed that the shift is 
along the (±l/8,7/8,0)-direction. The magnitude of the 
shift is proportional to the tilting angle and is assumed 
to be ~ 0.05a, 0.1a and 0.15a for the tilting angles of 
5, 10 and 15°, respectively, which are realisticjjalues for 
the compounds studied in the present workEO As for 
the Jahn- Teller distortion, it is assumed that the longest 
bond is by 0.1a longer than the shortest bond which is 
reasonable fat-LaMnOa and is relatively large for LaVOa 
and YTiOs.M 



III. RESULTS AND DISCUSSION 



A. LaMnOa 



system, in which one of the Cg or- 



In the high-spin d^ 
bitals is occupied at each site, the A-type antiferromag- 
netic (AFM) states with 'ix^ - r'^ / iy'^ 



dering with considerable mixture of 3z 



type orbital or- 
'^apipredicted 
to be stable by theoretical calculationsEl'Elil3 and are 
studied by x-ray and neutron diffraction measurements. □ 
Here, the z-direction is along the c-axis. The amount of 
the 3z^ — r^ component decreases with the JT distortion.^ 
Different ways of stacking the orbitals along the c-axis 
give two types of orbital ordering: the one compatible 
with the d-type JT distortion and the other with the a- 
type JT distortion. These two types of orbital ordering 
are illustrated in Fig. ||. While, in the orbital ordering 
of the a-type, the sites 1, 2, 3, and 4 are occupied by 



3y^ 



3x^ 
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sites 1, 2, 3, and 4 are occupied by iy 
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orbitals in the orbital ordering 



compatible with the d-type JT distortion. 

In Fig. |[ the energy difference between the orbital 
ordered states compatible with the d-type and a-type JT 
distortions is plotted as a function of the tilting angle of 
the octahedra, i.e., the magnitude of the GdFeOs-type 
distortion. A, C/, and (tWc) are 4.0, 5.5, and -1.8 eV, re- 
spectively, for LaMnOa.El Without the JT distortion and 
the shift of the ^-site ion, the two states are degener- 
ate within the accuracy of the present calculation (± 1 
meV/formula unit cell). This degeneracy is lifted when 



the JT distortion is included. With the JT distortion and 
the shift of the ^-site ion, the orbital ordered state with 
the d-type JT distortion is lower in energy than that with 
the a-type JT distortion. If we tentatively switch off the 
shift of the A-site ion and include only the JT distortion, 
the orbital ordered state with the a-type JT distortion 
becomes slightly lower than that with the d-type JT dis- 
tortion as shown in Fig. y. Therefore, one can conclude 
that the shift of the A-site ion driven by the GdFeOa-type 
distortion is essential to stabilize the orbital ordered state 
with the d-type Jahn- Teller distortion.EJ 

The qualitative explanation of this behavior is as fol- 
lows. In the d-type JT distortion, the four oxygen ions 
nearest to the A-site ion [shaded in Fig. n^ (a) and (b)] 
shift approximately in the same direction and, conse- 
quently, the system can gain the hybridization energy 
between the A-site and oxygen ions effectively. On the 
other hand, in the a-type JT distortion, the two of the 
four oxygen ions move in the other direction and the 
energy gain due to the hybridization is small compared 
to the d-type JT distortion. Another possible picture is 
that, in the d-type JT distortion, these four oxygen ions 
can push the A-site ion in the same direction since the 
JT distortion along the c-axis is in phase. On the other 
hand, in the case of the a-type JT distortion, the two 
oxygen ions in the upper plane push the A-site ion in the 
other direction than the two in the lower plane as shown 
in Fig. ^ Therefore, the stronger is the GdFeOs-type 
distortion, the more does it stabilize the d-type Jahn- 
Teller distortion and corresponding orbital ordering. 

In the charge-ordered state of Pro.sCao.sMnOs, the 
Mn'^+ and Mn^+ sites are arranged like a checkerboard 
within the c-plane_and the same arrangement is stacked 
along the c-axisJlj The Mn'^+ sites are accompanied 
by the JT distortion and the elongated axes are par- 
allel along the c-axis just like the d-type JT distor- 
tion. Since the A-sites are occupied by Pr and Ca 
ions in Pro.sCao.sMnOs, we cannot simply apply the 
present model calculation to it. However, it is reason- 
able to speculate that the stacking along the c-axis in 
Pro.sCao.sMnOa is also determined by the interaction be- 
tween the JT distortion and the shift of the A-siie ion in 
the same way as in LaMnOs. 



B. YVO3 

In the (P system, the C-type AFM state in which the 
sites 1, 2, 3, and 4 are occupied by xy and yz, xy and 
zx, xy and zx, and xy and yz orbitals and the G-type 
AFM state in which the sites 1, 2, 3, and 4 are occupied 
by xy and yz, xy and zx, xy and yz, and xy and zx 
orbitals are competing. While the C-type AFM state 
is favored by the orbital ordering which is compatible 
with the a-type JT distortion, the G-type AFM state 
is favored by the orbital ordering of the d-type. The 
relative energy of the G-type AFM state with the d-type 



JT distortion to the C-type AFM state with the a-type 
JT distortion, Ed — Ea, is plotted as a function of the 
tilting angle in Fig. ||. A, U, and {pda) are. 6.0, 4.5, and 
-2.2 eV, respectively, for LaVOs and YVO3.I Without the 
GdPeOa-type distortion, the C-type AFM state with the 
a-type JT distortion is lower in energy than the G-type 
AFM state, indicating that the energy gain due to the 
orbital ordering is larger in the C-type AFM state than 
in the G-type AFM state. The energy difference becomes 
smaller with the tilting or the GdFeOa-type distortion. 
Finally, with the tilting of 15°, the G-type AFM state 
with the d-type JT distortion becomes lower in energy 
than the C-type AFM state. The present calculation is 
in good agreement with the experimental result that the 
less distorted LaVOs is C-type AFM below 140 KQ aijud 
the more distorted YVO3 is G-type AFM below 77 K.B 
This situation is illustrated in Fig. g. When the 
GdFeOs-type distortion is large, the interaction between 
the d-type JT distortion and the shift of the A-site ion, 
namely, the energy gain due to A-0 covalency becomes 
dominant just like in LaMnOa and, consequently, the 
G-type AFM with the d-type JT distortion is favored. 
On the other hand, when the GdFeOa-type distortion is 
small, the energy gain due to orbital ordering becomes 
dominant and the a-type JT distortion and the orbital 
ordering compatible with it are realized. The JT distor- 
tion may be suppressed if the system is located near the 
crossing point where the a-type and d-typc JT distortions 
are almost degenerate. An interesting experimental re- 
sult related to this point is thai. YVO3 becomes C-type 
AFM between 77 K and 118 K.Q The present model cal- 
culation suggests that, if the JT distortion is switched off, 
the C-tyne AFM state is favored because of the orbital 
ordering.H Therefore, YVO3 is expected to be close to the 
crossing point and become C-type AFM when the d-type 
JT distortion is suppressed at elevated temperature. In 
this sense, between 77 K and 118 K, YVO3 may be an 
ideal orbitally ordered system without JT distortion. 



7.0, 4.0, and -2.2 eV, respectively, for YTiOs.B With the 
tilting of 15°, the orbital ordered state of the d-type is 
lower in energy, in agreement with the experimental re- 
sult. In this state, the sites 1, 2, 3, and 4 are occupied by 
ciyz + C2xy, cizx + C2xy, ciyz — C2xy, and cizx — C2xy or- 
bitals (ci ~ 0.8 and C2 ~ 0.6). However, experimentally, 
less distorted LaTiOa has no cuivery small JT distortion 
and has a G-type AFM stateJl3 The present calculation 
cannot explain why the G-type AFM state can be stable 
compared to the FM state in LaTiOa. 



IV. CONCLUSION 

In conclusion, we have studied the relationship between 
orbital ordering and the JT and GdFeOa-type lattice 
distortions. It has been found that the covalency be- 
tween the ^-site cations and oxygen makes the d-type 
JT distortion (same orbitals along the c-direction) lower 
in energy than the a-type JT distortion (alternating or- 
bitals along the c-direction) in the presence of the large 
GdFeOs-type distortion. As a result, the orbital ordered 
states compatible with the d-type JT distortion are fa- 
vored in LaMnOs, YVO3, and YTiOa which have the 
relatively large GdFeOa-type distortion. On the other 
hand, in less distorted LaVOs, the orbital ordering com- 
patible with the a-type JT distortion is favored because 
of the pure superexchange effect. 
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C. YTiOa 



For the d^ system, the ferromagnetic (FM) states with 
orbital ordering are favored in the model HF calcula- 
tions. There are two possible orbital orderings compati- 
ble with the a-type and d-type JT distortions. The model 
HF calculation without the covalency between the j4-site 
cation and the oxygen ion predicted that the orbital or- 
dering of the a-type is lower in energy. However, the 
recent neutron diffraction measurement by Akimitsu et 
al. have shown that the orbital ordering compatible with 
the d-type JT distortion is realized in the FM insulator 
YTiOa la YTiOa has the considerable GdFeOa-type dis- 
tortion and the tilting angle is expected to be larger than 
15°. In Fig. M, the relative energy of the FM and orbital 
ordered state of the d-type to that of the a-type is plot- 
ted as a function of the tilting. A, U, and (pda) are 
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FIG. 2. Two types of orbital ordering for LaMnOs com- 
patible with the a-type and d-type JT distortions. 
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FIG. 1. Schematic drawings for two types of Jahn- Teller 
distortions (a) and for GdFeOs-type distortion (b). The ar- 
rows indicate shifts of the oxygen and A-site ions. In the 
GdFeOa-type distortion, the six oxygen ions nearest to the 
A-site ion are shown. 



FIG. 3. Energies per formula unit cell of the orbital or- 
dered state compatible with the d-type JT distortion rela- 
tive to that with the a-type JT distortion for LaMnOa as 
a function of the tilting angle, i.e., the magnitude of the 
CdFeOa-type distortion, closed circles: with the JT distor- 
tions and with the shift of the A-site cation, open diamonds: 
without the JT distortions and without the shift of the A-site 
cation, closed triangles: with the JT distortions and without 
the shift of the A-site cation. 
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FIG. 4. A rough sketch of the shifts of the oxygen ions 
and the forces to the A-site ion for the a-type and d-type JT 
distortions. The oxygen shifts due to GdFeOa-type distortion 
are exaggerated. Only the shifts and forces relevant for the 
present discussion are shown. In the a-type JT distortion, the 
oxygen ions in the upper and lower planes push the A-site ion 
in different directions. In the d-type JT distortion, the oxygen 
ions push the A-cite ion in the same direction. 
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FIG. 5. Energies per formula unit cell of the orbital or- 
dered state compatible with the d-type JT distortion (closed 
circles) relative to that with the a-type JT distortion (open 
circles) for YVO3 and LaVOa as a function of the tilting. The 
shift of the A-site cation is included. 



FIG. 7. Energies per formula unit cell of the orbital or- 
dered state compatible with the d-type JT distortion (closed 
circles) relative to that with the a-type JT distortion (open 
circles) for YTiOa as a function of the tilting. The shift of 
the j4-site cation is included. 
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FIG. 6. Energy scheme of the orbital ordered states with 
the d-type and a-type JT distortions as a function of tilt- 
ing or the GdFeOa-type distortion. While, for the small 
GdFeOa-type distortion, the energy gain due to the orbital 
ordering is dominant, the energy gain due to the A-0 cova- 
lency becomes relevant for the large GdFeOa-type distortion. 



